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ABSTRACT 

We perform an analysis of ~ 80000 photometric measurements for the following 10 
stars hosting transiting planets: WASP-2, -4, -5, -52, Kelt-1, CoRoT-2, XO-2, TrES-1, 
HD 189733, GJ 436. Our analysis includes mainly transit lightcurves from the Exo¬ 
planet Transit Database, public photometry from the literature, and some proprietary 
photometry privately supplied by other authors. Half of these lightcurves were ob¬ 
tained by amateurs. From this photometry we derive 306 transit timing measurements, 
as well as improved planetary transit parameters. 

Additionally, for 6 of these 10 stars we present a set of radial velocity measurements 
obtained from the spectra stored in the HARPS, HARPS-N, and SOPHIE archives 
using the HARPS-TERRA pipeline. 

Our analysis of these TTV and RV data did not reveal significant hints of addi¬ 
tional orbiting bodies in almost all of the cases. In the WASP-4 case, we found hints 
of marginally significant TTV signals having amplitude 10 — 20 sec, although their 
parameters are model-dependent and uncertain, while radial velocities did not reveal 
statistically significant Doppler signals. 

Key words: planetary systems - techniques: photometric - techniques: radial veloc¬ 
ities - methods: data analysis - methods: statistical - surveys 
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1 INTRODUCTION 

The first extrasolar plan et, orbiting a sola r -type star 51 Pe- 
gasi, was discovered by iMavor fc Queloa 1 1995 1, based on 
the precision Doppler measurements of the ELODIE spec¬ 
trograph. After that, the number of the detected exoplanets 
grew continuously, exceeding 1000 so far. In fact, 20 years 
ago a new rapidly-growing domain of fundamental science 
was created, devoted to the exoplanet research. 

Currently, most of the known exoplanetary candidates 
were detected by one of two major techniques: the radial 
velocity (RV) method or transit method that increased its 
output in recent time thanks to the launch of specialized 
spacecraft CoRoT (ESA) and Kepler (NASA). 

Unfortunately, both the Doppler and transit exoplanet 
detection methods require sophisticated and expensive in¬ 
strumentation and remain practically inaccessible to the ma¬ 
jority of the community. The RV method requires the use of 
extremely stable spectrographs. There is only about a dozen 
of such instruments in the world that are capable of detect¬ 
ing an exoplanet. The transit detection is also difficult. It 
requires precise alignment of the planetary orbit with the 
observer, and the transit event is rather short in time, al¬ 
though periodic. To detect a transiting planet, we have to 
observe lots of stars, and we necessarily deal with large num¬ 
bers of null detections. This requires the use of specialized 
ground-based robotic telescopes or telescope networks, or 
space observatories like the above-mentioned CoRoT and 
Kepler. Organizing such a campaign is a difficult task even 
for professional scientific teams. 

In terms of the photometric accuracy and quality, 
ground based telescopes are no match to spacecrafts like 
Kepler. But all space projects have a common disadvan¬ 
tage: they are severely limited in time, while the sensitiv¬ 
ity to weak planetary signals in the data degrades quickly 
when the planetary orbital period exceeds the observational 
time base. This condition means that Kepler data cannot 
reliably detect long-period planets, analogous to the giant 
planets of Solar System. In general, the responsibility for 
follow-up observations always returns to ground-based ob¬ 
servatories. Ground-based observations of planetary transits 
are much less demanding than precision radial velocity mea¬ 
surements, and of course much cheaper than projects like 
Kepler. In fact, such observations are possible with com¬ 
mercially available equipment typically used by a significant 
community of “amateurs”. Thus, the exoplanetary hunt can 
potentially become “citizen science” that can trigger a qual¬ 
itative leap in the field. 

However, amateurs are definitely not equipped to un¬ 
dertake classic transit surveys like e.g. SuperWASP or oth¬ 
ers. But nonetheless they can provide a useful scientific con¬ 
tribution b^jneans_of_thetrm^t timirig_variation (TTV) 
method llHolman fe Murravll2005l : lAgol et al.ll2005h . In this 
approach, there is a list of well defined targets that are con¬ 
tinuously monitored. Each target is a known host of a tran¬ 
siting planet, and its transits are regularly observed. If more 
planets orbit the host, they should induce perturbational ef¬ 
fects on the motion of the transiter, causing observable de¬ 
lays to its transits, in comparison with a strictly periodic 
ephemeris. The TTV method allows the observation time 
to be used more efficiently, because we know what targets 
should be observed, and when. 


This work represents an attempt to determine the prac¬ 
tical efficiency of such an approach, based on the pho¬ 
tometry data, taken mainly from the Exoplanet Tran¬ 
sit Database (ETD) of the Czech Astronomical Union, 
http://var2.astro.cz/ETD/. 

Currently, about 30 observatories are regularly con¬ 
tributing to the database, including amateur as well as pro¬ 
fessional ones. This network offers telescopes of different 
apertures - from 20 cm to 2.6 m. Their locations are shown 
in Fig.Q] We must note that in terms of technical character¬ 
istics of the telescopes, there is no sharp boundary between 
the amateur and professional equipment, and the quality of 
the observations is often determined by the local astrocli- 
mate, which can have even more important effect than the 
telescope size. The observational programme involves cur¬ 
rently ~ 20 transiter hosts that are more or less regularly 
observed. For some of the stars, several years of observations 
are already available. However, the transit fitting algorithm 
used by E TD is criticized for its simplicity and imperfec¬ 
tions (e.g. IPetrucci et al.l [20131 '). In this work we present a 
new data reduction pipeline that handles subtle photomet¬ 
ric effects (like the red noise) more accurately. This pipeline 
was developed to deal with the photometry of a relatively 
poor or moderate quality, which is typical for the amateur 
data in ETD. 

We do not limit ourselves to only amateur observations 
or only ETD. We also use photometric data published in the 
literature. Moreover, for our targets we revealed a moder¬ 
ate amount of the spectra stored in the HARPS, HARPS-N, 
and SOPHIE archives, and we derive RV data from them to 
provide and independent “calibration” of our TTV results. 
However, in this work we have no goal to perform a fully 
self consistent TTV+RV analysis. Instead, we aim to char¬ 
acterize the accuracy and reliability of the TTV data that 
we can derive from just the photometry. 

The structure of the paper is as follows. In Sect. [2] we 
provide a detailed description of all the data that we in¬ 
clude in our analysis. In Sect. [ 3 ] we introduce the algorithms 
used to process the photometric data. In Sect. [4] we present 
the TTV data derived from the photometry and describe 
the results of their periodogram analysis. In Sect. [5] we give 
the remaining fitted parameters of the planets considered in 
the work. In Sect. [S] we describe the RV data obtained for 
some of our targets on the base of the spectra found in the 
HARPS/HARPS-N/SOPHIE archives. In Sect.[7]we discuss 
in detail the case of WASP-4, for which we detected possible 
hints of a weak TTV signal. 


2 THE SOURCE DATA 

Our primary source of the photometric data is the Exoplanet 
Transit Database (ETD). We extracted a set of the best 
transit lightcurves from ETD for 10 selected stars. Also, we 
added in our analysis several public transit lightcurves that 
we found in the literature. Thus, our photometric data cover 
10 targets with 306 transit lightcurves, and about 80000 
individual photometric measurements. The public data were 
presented in the works listed in Table [T] and most are stored 
in the Vizier database. Among these 306 transit lightcurves, 
161 (or roughly half) were obtained by amateurs, while 65 
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Figure 1. World distribution of the observatories regularly contributing to the ETD database. 


Table 1. Sources of the photometric data (except for ETD). 


Table 2. Explanation of the RV data file 


target 


references 


WASP-2 

WASP-4 


WASP-5 

WASP-52 

Kelt-1 

CoRoT-2 


Wilson et al. 

2008; Gillon et al.l |2009 

Winn et al. 

200E& 

Southworth et al. 2009tj 

Sanchis-Ojeda et al. 

201ll;lNikolov et al. 2012 


Petruccietal. |2013| 


Southworth et ai~li2009al 


No 


Siyerdetal. 2012 


Gillon et al. 201C 


TrES-1 

Winn et al. 

2007b] 

HD 189733 

Bakos et al. 

,2006; 

IWinn et al. 

200H 


Pont et al. 2007 

GJ 436 

Bean et al. 

m 

o 

o 

CM 

Shoorer et a] 

2009 


Stevenson et al.ll2012l 


are from professional observatories contributed to ETD, and 
80 were published in the listed literature. 


We need to note that the data from 


(Sanchis-Oieda et al. 

20 lib also contain reprocessed 

photometry from 

Winn et al. 2009), so the original 


Iwinn et al] (2009) data were not actually included in our 
analysis. Als o , we e ventu ally decided to drop the HST 
data by iBean et al.l (120081 ). because they all appeared to 
cover only small parts of a transit and could not produce 
good precisio n in the derive d mid- times. The photometry 
presented by IPetrucci et ali (|2013l ) is not public, but it 
was kindly released to us by the au thors. The WASP-4 
data from (ISouthworth et all l2009bl ) appeared not very 
reliable due to the clock errors of the Danish telescope 
dNikolov et al.1120121 : IPetrucci et al.ll20 13l). a nd we believe 
that the data from (ISouthworth et al.ll2009al . l20ld ) should 


column explanation 

1 observation time (BJD) 

2, 3 derived radial velocity and its uncertainty 
4 a standardized name of the input RV file that includes 
the target name and the name of the spectrograph 


be treated with the same care. However this photometry 
is rather accurate in itself, and thus it is still useful in 
constraining all transit parameters except for mid-times 
(e.g. by adding a fittable time offset to these light curves 
model relative to the other lightcurves). All timings in the 
photometric data were transformed to the BJDtdb time 
stamps by means of the public IDL software developed by 
lEastman et al.1 (]2010l ). 

Additionally, we used the precision radial velocity (RV) 
obtained from the spectra of the HARPS archive, avail¬ 
able for the following targets from our photometry sample: 
WASP-2, -4, -5, HD 189733, Corot-2. These spectra were 
processed with the advanced HARP S-TER RA pipeline that 
offers an improved RV accuracy (|Anglada-Escud6 fc Butler! 
SH). The quality of these RV data, and their number per a 
target, appeared not very high. The meaning of the columns 
in the attached RV data file is given in Table [2] 

For GJ 436 we found a large amount of high- 
quality HARPS RV data, similar to the one considered by 
lLanotte et al.l (120141 ). as well as a larg e Keck RV time series 
recently published by iKnutson et al.l J2014I ). We acknowl¬ 
edge that cases like this deserve to be investigated detailedly 
in a separate work. In this paper we only present our TTV 
and TERRA RV data for GJ436, processed in the common 
simplistic way as for the other stars. 

In this work we do not perform the joint transits+RV 
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fits, because it appeared that to perform such an analysis 
at a desirable level of quality, we must provide a solution 
to several non-trivial issues that fall outside the scope of 
this paper. These issues include, e.g., the treatment of the 
correlational structure of the RV noise, and treatment of 
the RV points affected by the Rossiter-McLaughlin effect. 
In this work we used the RV data mainly as an independent 
source of information. 


3 METHODS OF THE TRANSIT 
LIGHTCURVES ANALYSIS 


3.1 Details of the transit model and its 
parametrization 


Let us first adopt the approximation that the transiting 
planet moves along a straight line with constant velocity, 
thus neglecting the orbital curvature of its actual trajectory 
during the transit. The model of the linear motion easily pre¬ 
dicts the separation between the centers of planet and star 
disks, 5, as a function of time. We have 4 kinematic charac¬ 
teristics of the transit that must be fitted: (i) the mid-time 
of the transit t c , (ii) the duration of the transit td, defined 
as the time spent between the first and fourth contacts and 
(iii) the impact parameter b, measuring the smallest pro¬ 
jected separation <5, and (iv) the projected planet radius r 
that simultaneously determines the transit depth and the ge¬ 
ometry of the ingress/egress phases. Given these parameters 
and easy geometric constructions, the projected separation 
S can be expressed as: 

5(t) = \/b 2 + [(1 + r) 2 — b 2 ] t 2 , r = 2-—— , |6| ^ l+r.(l) 

td 


This assumes that the star radius is unity. To simplify the 
formula for 5 and simultaneously get rid of the non-trivial 
definition domain for b, our algorithm adopts internally a 
replacement of the impact parameter b: 




:(1 +r). 


( 2 ) 


The trajectory curvature also depends on the orbital 
eccentricity, but usually this eccentricity is difficult to de¬ 
termine from transit observations reliably, and we have no 
other option except to assume that it is zero. A zero eccen¬ 
tricity is a good prior assumption for most of these short- 
period planets. In the cases when radial velocity data are 
also available, the accurate eccentricity information can be 
obtained from the joint transit+RV fits, but we do not per¬ 
form fits of such type in this work. 

After the function 5(t) is calculated, we approximate 
the relative flux red uction with the use of the stel lar limb 
darkening model by lAbubekerov fc Gostevl (l2013l j. These 
authors provided theoretical formulae for various types of 
the limb darkening effect, as well as a software library writ¬ 
ten in C. The library provides subroutines to compute the 
observed light flux reduction A L as a function of the eclips¬ 
ing planet radius r and of the projected separation 5 (as¬ 
suming that the star radius is unit), as well as of the limb 
darkening coefficients that depend on the selected model. 
Besides, this library provides partial derivatives of AL with 
regard to its arguments. These derivatives are necessary to 
compute the gradient of the likelihood function, which is 
also used by the transit fitter. 

Thus for a transit lightcurve we have 4 fittable param¬ 
eters of the planet or planetary orbit: the planet/star radii 
ratio r, the transit mid-time t c , the transit duration td, and 
the impact parameter replacer p. As the orbital period P 
is known for all our planets with a very good accuracy, we 
treat it as a fixed parameter in the formula 0. 

Additionally, there are two parameters determining the 
limb darkening model. In this work we use a quadratic two- 
term model of the stellar limb darkening with two coeffi¬ 
cients to be determined, A and B. The brightness of a point 
on the visible stellar disc, observed at a given separation 
from its center, p, is modelled as 

I(p) = 1 — A(1 — p) — B(1 — p) 2 = 

= 1 - A-2B + (A + 2B)p +Bp 2 , 

p = \/l -p 2 . (5) 


With the new parameter p we do not need to worry about its 
domain: any real value of p is physically meaningful. Thus 
we eliminate the danger that this parameter can walk to a 
forbidden domain during the fit. With p replacing b, we have 


5{t) = (1 + r) 


/ 


T 2 _|_ p2 

1 +P 2 ' 


( 3 ) 


However, this approximation of 8(t) might be inaccu¬ 
rate due to the curvature of the planet trajectory. In this 
work we include a correction to this formula, assuming that 
the planet moves along a circular orbit. In this approxima¬ 
tion, the projected distance can be expressed by the same 
formulae 0 and ©, substituting the following quantity t' 
instead of r: 


sin to: 
sin a ’ 



« 


1 , 


( 4 ) 


where P is the transiter’s orbital period. The auxiliary angle 
a reflects the curvature of the circular orbit. In fact, this 
correction induced only a well negligible effect on our TTV 
measurements, but nonetheless our results correspond to the 
model with this correction included. 


However, in practice the brightness model @ can easily 
turn non-physical, if the coefficients A and B are allowed to 
attain arbitrary values. This becomes a significant problem 
if our data are polluted by some systematic errors that are 
always difficult to foresee in advance. There are a couple of 
natural basic constraints that we place on the coefficients A 
and B to keep the model © physically reasonable: 

I{p)> 0, -y- = [2Bp — (A + 2B)\ — ^ 0, Vp € [0,1].(6) 

dp p 

If the second condition (the one on the derivative) is sat¬ 
isfied, 7(p) is a monotonic non-increasing function, so the 
condition 7(p) ^ 0 for p £ [0,1] is then equivalent to 
7(1) = 1 — A — B ^ 0. In the second condition, the ex¬ 
pression in parenthesis is a linear function of p, so to have 
it always non-positive for p £ [0, 1], it is necessary and suf¬ 
ficient to have it non-positive at the boundaries p = 0 (disk 
limb) and p = 1 (disk center). Finally, we have total of three 
elementary inequality constraints on A and B to satisfy: 

A + B^l, A + 2B>0, A> 0. 


( 7 ) 
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Note that the coefficient B is allowed to be negative here 
(but B > —1)- We do not put an extra condition B ^ 0. 
Negative values of B mean that the limb-darkening gradi¬ 
ent is diminishing (in absolute value) closer to the limb, 
although it never turns positive (any “limb brightening” is 
disallowed). 

We may satisfy 0 by means of making a smooth re¬ 
placement of the parameters A and B, such that the formu¬ 
lae of the replacement would disallow the conditions 0 to 
be broken. This can be reached, for example, by the follow¬ 
ing trigonometric replacement: 

A = sin 2 9 (1 — cos p), B = sin 2 9 cos p. (8) 

Naturally, whatever real values the new auxiliary parameters 
9 and p might attain, the resulting values of A and B always 
satisfy 0 . From the other side, each point {A, B) in the 
domain 0 maps to some pair (9 , p) with real values of the 
parameters: 

9 B 

sin 2 9 = A + B, cosi p =—. -—. (9) 

A + B v ’ 

Note that from © it follows that A + B ^ |B|, meaning 
that A + B is never negative. 

Therefore, treating the auxiliary angles 9 and p as pri¬ 
mary fittable parameters, we can satisfy the conditions 0 
automatically. The result of the fitting would correspond 
to either an internal point of the domain 0 , or to some 
point on its boundary, if the actual data suggest to move 
the solution to a non-physical domain due to e.g. their poor 
statistical quality or systematic errors. The boundary points 
correspond to certain special values of 9 and p that can be 
easily identified: 9 = ±7r/2 + 7rfc on the first boundary of 0 , 
p = 7T + 2nk on the second boundary, and p = 2nk at the 
third boundary. 


3.2 Details of the fitting procedure 

In our analysis we are interested in obtaining the TTV devi¬ 
ations with maximum accuracy achievable with the available 
photometric data. However, the ETD data, even after the 
pre-selection, often have only a moderate quality, performing 
a complete and independent fit for each transit lightcurve is 
not a good option. The accuracy of thus obtained transit pa¬ 
rameters would often be poorly constrained, and this would 
impact the accuracy of the fitted TTV data as well. More¬ 
over, the transit fits for individual lightcurves sometimes do 
not even converge to a reasonable result, e.g., due to incom¬ 
plete coverage of the transit or presence of significant curved 
trends. 

To overcome these issues, we adopt in this work the 
following approach. We perform a joint fit of all transit 
lightcurves, available for a given star, assuming that most of 
the transit parameters, except for the mid-times, are equal 
for different lightcurves. Such a “shared” transit parameter 
is still fittable, taking into account the constraint of its val¬ 
ues being equal between different lightcurves. The mid-times 
are fitted individually for each lightcurve, i.e. they remain 
unconstrained. 

Such an approach uses the full statistical power of all 
photometric data, available for a given star, to fit the shape 
of the transit curve, while still fitting the TTV offsets of 


individual transits separately from each other. We can note 
the following potential weaknesses of this approach: 

(i) It is not taken into account that the limb darkening 
coefficients depend on the photometry band, which are dif¬ 
ferent for different lightcurves. Therefore, the results of such 
fitting method would refer to some averaged limb darken¬ 
ing, possibly introducing minor modelling errors in individ¬ 
ual lightcurves. From the other side, the limb darkening ef¬ 
fect in the transit curve is always symmetric relative to the 
mid-time, implying that its impact on the derived mid-times 
should be small, if not negligible. In fact, we noticed that 
even fitting of a transit model without any limb darkening 
at all does not change the derived TTV data significantly 
(beyond the estimated parametric uncertainties). However, 
in the final version of our algorithm we decided to disentan¬ 
gle the limb darkening coefficients for the best lightcurves 
(those that have r.m.s. smaller than 10 per cent of the transit 
depth) and fit these coefficients independently. 

(ii) It is not taken into account that transit duration 
may also be subject to variations, like the transit mid-time. 
However, in practice, it appeared that the accuracy of the 
transit duration estimations was roughly an order of mag¬ 
nitude worse than those of the mid-times. Moreover, some 
lightcurves cover the transit event only partially, implying 
that its duration would remain almost unconstrained when 
fitted independently. We do not address transit duration 
variations (TDV) in this work, focusing our attention on 
the TTV. 

(iii) The mid-time estimates, obtained in such a manner, 
are not necessarily uncorrelated. These mid-times are corre¬ 
lated with the remaining transit parameters, which are now 
shared between the lightcurves. Through this effect, some 
cross-correlation between the estimated mid-times may ap¬ 
pear. This creates a risk of unexpected statistical effects in 
the derived TTV data, like e.g. the non-white noise. How¬ 
ever, the magnitude of these TTV correlations usually re¬ 
mains very small (maximum of a few per cent, and ~ 0.1 per 
cent in average), and no deviations from the white noise are 
seen in the periodograms of the TTV data. More significant 
(> 10 per cent) TTV correlations can appear for lightcurves 
that offer only a partial coverage of the transit, but such 
lightcurves represent only a minor fraction of our data. 

In addition to the planetary transit itself, our lightcurve 
models include a polynomial trend with fittable coefficients. 
Such a trend is necessary to take into account various drift¬ 
ing effects, e.g. the effect of airmass or other types of sys¬ 
tematic variations that appear frequently in our data. Each 
transit lightcurve has an individual fittable trend with a 
separate set of trend coefficients. After some experiment¬ 
ing with the data, we found that cubic trends represent a 
good compromise between the model adequacy and its para¬ 
metric complexity. We did not try to reduce this systematic 
photometric variation based on its correlation with the air- 
mass function: while for some lightcurves such a correlation 
looked clear, for others it was not obvious, indicating that 
other systematic effects were in the game. 

The lightcurves model is relatively complicated, while 
the photometry used in this work is not of a very good qual¬ 
ity. In practice we often faced difficulties causing the fitter 
to be trapped in the local maxima of the likelihood function, 
which was related to an unrealistic branch of the solutions. 
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For example, without special care, see below, we frequently 
obtained non-physical solutions corresponding to a grazing 
transit with r 1. Also, some lightcurves do not cover the 
complete transit, and in such cases the fitting of the transit 
mid-time is complicated by its strong correlation with the 
coefficients of the polynomial trend. To avoid such traps, we 
worked out the following sequence of auxiliary preliminary 
fits: 

(i) Perform a preliminary fit constraining the mid-times 
at a regular grid (with free scale and offset), implying all 
TTV residuals are zero by definition; fixing the impact pa¬ 
rameter p at an intermediary value of 1 (6 « l/\/2), and 
fixing the limb darkening coefficients at A = B = 0.25 (or 
9 = 7t/4 and y? = 7t/3). 

(ii) Refit after releasing the mid-times and impact param¬ 
eter, but still holding the limb darkening coefficients fixed. 

(iii) Refit after releasing the limb darkening coefficients 
(binding them across different lightcurves). 

(iv) Refit after full release of the limb darkening coeffi¬ 
cients for the best lightcurves (those with r.m.s. < 0.1 of 
the transit depth). 

After that, we also apply the red-noise detection and 
fitting procedure as described in the section below. 


3.3 Reduction of the red noise 

It is already known well that stellar photometric data 
usually include a correlated (“red”) noise component that 
has an importan t effect on the fitted transit parameter s 
llPont et al] 120061 : IWinn et al.l 120081 : ICarter fe Winnl l2009h . 
An easy technique to calculate the effect of the red noise on 
the fitting uncertainties was introduced in these works. Al¬ 
though these authors avoid making restrictive assumptions 
concerning the correlation structure of the red noise, their 
method still remains rather simplistic and it does not take 
into account important effects. They mainly focus on a more 
accurate determination of the uncertainties in the transit fit, 
which are typically underestimated without a proper treat¬ 
ment of the red noise. But the red noise may also induce 
biases in the fitted parameters themselves. This biasing ef¬ 
fect was already noted when processing Doppler data af¬ 
fected by re d noise in a generally simialar way (e.g. iBaluevI 
l201ll. l2013blh Also, there is not an obvious way to control 
the validity of the underlying assumptions in the traditional 
approach of the photometric red noise reduction. This ap¬ 
proach does not offer a complete noise model that could be 
verifi ed for accur acy. Additionally, the original method by 
IPont et alj (l2006li relies on the out-of-transit photometric 
data, which are very limited in our case. 

In this work we apply the approach based on a para¬ 
metric modelling of the noise covariance matrix. Th is is an 
adap t ation of the red-noise fitting technique from (IBaluevI 
l201ll . l2013bl . 120151 1 that was developed to handle exoplan¬ 
etary radial-velocity data, in which the noise correlations 
may appear e.g. due to the stellar activity. In this approach 
we approximate the red noise by a stationary Gaussian ran¬ 
dom process with a covariance function of a given functional 
form. The covariance coefficient between two arbitrary pho¬ 
tometric observations Xi and Xj, acquired at the times ti 


and tj, is modelled as 

I/y — Gov (Xj , Xj ) — (7 , ^(pv.'ht )Gy T V/j',red (pred , v) , 

Vij, red =.Pred Rij{r), Rij = p ^ ^ ^ , (10) 

with p w ht, Pred, and r being free fittable parameters. Here, 
the first noise term, er? wht , represents the white fraction 
of the noise, detailed below. The red noise is given by the 
second term, Vij, which also depends on the function p(t) 
that represents an adopted shape of the correlation func¬ 
tion. We use mainly the exponential correlation function 
p(t) = exp(—|f|). This not a unique choice, but in practice 
this model of the red noise usually appears adequate. 

The white-noise term erf wht in (HOD requires a separate 
discussion. We believe that physically the so-called additive 
model would be more suitable here. In this model, the value 
erf wht is determined as a sum of the stated instrumental vari¬ 
ance and of an unknown fittable “jitter” variance, generated 
by e.g. Earth atmospheric instability, or by unassessed in¬ 
strumental instability, or by short-term (minutes to hours) 
intrinsic variations of the stellar flux, whenever these vari¬ 
ations can be treated as a white noise. Although, we must 
note that in practice the instrumental uncertainties in our 
data do not look very trustable, and often they are just 
omitted, forcing us to assume that the measurements have 
just equal uncertainties. Therefore, in any case we should 
not expect that our white noise can be accurately modelled 
from the physical point of view. In these circumstances, the 
noise model should be chosen mainly on the basis of math¬ 
ematical simplicity or usefulness, relying on only minimal 
or no knowledge of the underlying physics. We dec ided to 
use in our work a regularized noise model defined in (IBaluevI 
l2015l l. In the majority of practical cases, this model should 
be equivalent to the additive noise model. This regularized 
model proved rather resistant with respect to various pitfalls 
appearing during the fitting procedure. 

The fitting of the compound model, involving simulta¬ 
neously the models of the transit curve and of the photo¬ 
metric noise, is done by mean s of t h e max imum-likelihood 
approach with details given in (tBaluevil2015l i. That work was 
devoted primarily to the analysis of the radial velocity data, 
but mathematically the methods that we are using here are 
identical. 

However, after an attempt to apply this technique in 
practice literally, we faced the problem that the red noise 
could not be fitted in many of our lightcurves. Only 1/3 to 
1/2 of the lightcurves allowed for a reliable red noise fit, 
while in the other cases, the red noise was either not de¬ 
tectable, or its estimated magnitude became pretty small 
(in comparison with the uncertainty). Such cases lead to 
model degeneracies and decrease the reliability of the entire 
fit for a given star. Thus, we decided to add the red noise 
term to the photometric model only in the cases when it was 
justified. 

We added to our analysis pipeline a set of auxiliary fits 
in order to identify the lightcurves, in which the red noise 
could be modelled more or less reliably. Namely, we tried to 
perform a series of test fits, latterly adding a red noise term 
to the model of each transit lightcurve. If the red noise could 
not be estimated reliably, more accurately if the relative 
uncertainty of its estimated magnitude exceeded 2/3, the 
relevant red noise term was removed from the model and this 
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lightcurve was further treated as having purely white noise. 
Otherwise, we proceeded to the next test fit preserving this 
red noise term in the model. In the end, after all individual 
lightcurves were processed in such a way, we performed one 
more check of each red noise term (because some of them 
could turn insignificant after the other terms being added) 
and the final fit was made. 

The quality of such an approach to the red noise reduc¬ 
tion is examined in Fig. [3] The best way to observe a “non¬ 
white” noise is to look at its power spectrum: the white noise 
should have constant power in average, while the red noise 
should demonstrate a systematic uprise to long periods. In 
our case we com pute the so-called residual periodograms 
fsee iBaluevI 120151 1 for each individual lightcurve, and then 
compute their average per each star involved in the analy¬ 
sis. These averaged periodograms are plotted in Fig. [2] For 
each star, we sequentially try three base photometric mod¬ 
els: (i) just transits without any limb darkening plus the 
white noise; (ii) same as model (i) plus cubic trends and the 
effect of the quadratic limb darkening; (iii) same as model 
(ii) plus the red noise term. 

We can clearly see from Fig. [2] that the long term pho¬ 
tometric variations are not reduced to just trends. Although 
the trends themselves are also important, after their removal 
the power spectra still contain an additional power in the 
period range of > 10 min. The trends can only affect the 
period domain of > 200 — 300 min (this corresponds to the 
typical time span of the lightcurves). In the most cases, ap¬ 
plication of our red-noise model remarkably suppresses the 
remaining excessive power, at least by a factor of 2. We must 
acknowledge that the reduction of this red noise looks rather 
difficult and still not entirely complete. 

Some statistical information on the derived red noise 
characteristics now follows. About 60 per cent of lightcurves 
did not contain any detectable red noise at all. For the re¬ 
maining 40 per cent of lightcurves, the red noise caused an 
increase of the TTV uncertainty by a factor of 1.3 on av¬ 
erage (a median value). This is not very large, although a 
minor fraction of the mid-times demonstrated an increase in 
the uncertainty up to a factor of 2 — 3. Concerning the bias, 
appearing in the derived mid-times themselves (rather than 
uncertainties), its average value was only about 1/6 fraction 
of the uncertainty, and only a few points demonstrated a 
shift above la. 

Summarizing, in the majority of the cases, either there 
was not any detectable red noise, or the change in the 
TTV data was small, even if the red noise term was de¬ 
tectable. Nonetheless, it is important that some individual 
TTV points are affected much more. 

3.4 Testing the Gaussianity of the noise and 
clearing away the outliers 

In the data-analysis methods described above, we largely re¬ 
lied on the assumption that the photometric noise is Gaus¬ 
sian. This assumption needs some verification, which we per¬ 
formed in the following manner. After fitting all the data, 
we computed the residuals and normalized them by their 
relevant modelled uncertainty (square root of the sum of 
the estimated white and red noise variances). If the noise 
was Gaussian, the distribution of these normalized residuals 
should be close to the standard Gaussian, and vice versa. 


We consider the entire set of ~ 80000 normalized residuals 
for all stars involved in our analysis, and plot the relevant 
distribution in Fig. [3] 

As we can see, the empirical distribution is indeed al¬ 
most Gaussian, except for the very tails (> 4-sigma devia¬ 
tion). The empirical distribution has heavier tails. However, 
by removing data in the distribution tails we can reach a 
good agreement with the Gaussian approximation. In fact, 
these extreme points represent outliers and should be re¬ 
moved in any case. The final results below correspond to 
the data with these 24 outliers removed. We also notice that 
the minor systematic deviations still remai ning in the distri¬ 
bution tails are mainly due to iBakos et alj (l2006lf data for 
HD 189733. Visual investigation reveals that some of these 
lightcurves contain limited segments with une xpected ly in¬ 
creased photometric scatter. Removal of the Bakos et ali 
(I2006T) data makes the agreement of the residuals distribu¬ 
tion with the Gaussian one very good. 


4 DERIVED TTV DATA AND THEIR 

ANALYSIS 

Using the fitting techniqe described above, we computed the 
mid-times for each observed transit, and placed these TTV 
measurements and other accompanying data in the online 
supplement to the article. The supplement represents a file 
containing a text table. The meaning of the columns in this 
Hie is described in Table [3] This file only contains the pa¬ 
rameters that are attached to individual lightcurves. The 
common fit parameters that are shared between different 
lightcurves are given separately in Sect. [5] below. 

The TTV residuals for these data are plotted in Fig. [3] 
In Fig. [5] we verify the normality of the distribution of the 
TTV noise present in these derived data in the way simi¬ 
lar to Sect. nm In the TTV residuals we find no obvious 
outliers or deviations from the Gaussian distribution. The 
apparent excess in the distribution tails that can be seen in 
the left panel of Fig. [5] is uncertain due to a relatively small 
number of TTV data points. Nonetheless, we identified 9 
possible “candidate outliers”, and tried t o remove them to¬ 
gether wit h 7 T TV points referring to the lSouthworth et al.l 
(l2009alfbl. l2010l) data. The TTV analysis discussed below 
was performed for the full TTV data set as well as for the 
reduced one. 

First, we made an effort to detect possible long-term 
curvature in our TTV. This was done by means of modelling 
the mid-times by a quadratic function (instead of only a 
linear function that would correspond to strictly periodic 
transits). None of the TTV time series demonstrated a long¬ 
term quadratic variation with at least 2-sigma significance 
(the most suspicious cases had 1.5cr to 1.7<r). 

Second, we undertook a periodorgam analysis of the 
TTV residuals, in order to reveal possible periodic pertur¬ 
bations from any additional unknown planets orbiting these 
stars. We again applied “residual periodograms” similar to 
the ones that were used to plot Fig. [3] above. Now our null 
model included a quadratic trend that is related to the pe¬ 
riod of the main (transiting) planet and white noise with 
a fittable additive “jitter”. The alternative model also in¬ 
cluded a sinusoidal signal at a probe frequency. None of 
these periodograms demonstrated any signs of a non-white 
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Figure 2. Examining the effect of the red noise in the lightcurve photometry for a few of stars belonging to our sample. Please note 
that plots in the first column have different ordinate scale than those in the other two. See Sect. [3?3l for a detailed discussion. 
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raw data (78153 points) 24 outliers removed 




-10 -5 0 5 10 -10 -5 0 5 10 

normalized residual normalized residual 

Figure 3. Testing the normality of the entire ~ 80000 photometric data used in the work. Abscissa shows the value of the best fitting 
residual, normalized by its total modelled variance (for white+red noise). The ordinate contains the normal quantiles of the empirical 
cumulative distribution of these normalized residuals. Perfectly Gaussian residuals should all lie close to the main diagonal that represents 
a standard normal distribution, while a deviation from the diagonal indicates a non-Gaussian noise. 


noise, and in fact almost all appeared consistent with the 
pure noise model, although the cases revealing significant 
excessive TTV “jitter” were rather frequent. No significant 
TTV periodicities were found in the data for any of the stars 
involved in the analysis, except for the WASP-4 case which 
is discussed below. The periodog ram sign i ficance leve ls were 
estimated using the approach of lBaluevl 120081 . [20091 ). 

We believe that this result indicates a considerable im¬ 
provement in the statistical quality of the TTV data and/or 
methods of the statistical analysis, because for the original 
TTV data from ETD we frequently detected spurious pe¬ 
riodicities possibly appearing due to imperfections of the 
transit fitting algorithm. 


5 IMPROVED TRANSIT CURVE 
PARAMETERS 

In addition to the TTV data, we also provide the re¬ 
maining best fitting parameters of the exoplanetary tran¬ 
sit lightcurves. These are the transit parameters that were 
shared between different lightcurves, and they are given in 
Table [4] In addition to the previously mentioned transit pa¬ 
rameters r, td, b, and the limb darkening parameters A and 
B, this table also contains the following data: number of 
the transit lightcurves used in the analysis, number of the 
lightcurves in which a red noise term was robustly detected, 
the x 2 (weighted r.m.s.) of the TTV residuals (relatively to 
the best linear fit of the derived mid-times), and the refined 
parameters of the transit emphemeris: transiter period P 
and a reference mid-time t() ef . The last two quantities were 
obtained by requiring the mid-times to obey a strict linear 
relation with the transit count and fitting the coefficients 


of this relation. The value of t/ ef is not unique due to the 
periodic nature of the transits. We chose such a reference 
transit, for which the uncertainty of t c would be minimum. 
This condition simultaneously implies that the correlation 
between the estimated P and t() ef is zero. Note that it is not 
required to have actual transit observation at this t/ ef . 

The Xttv values in Table [4] always exceed unity, so 
the fit uncertainties might be moderately underestimated. 
In some cases this higher than expected scatter could be 
partly explained by small number of transits involved in the 
analysis, and in some part it might be due to additional un¬ 
seen planets that may induce complicated TTV signals and 
are difficult to detect. Another explanation is incomplete re¬ 
duction of the red noise and the effect of nonlinearity of the 
transit model. The latter effect can introduce biases both 
in the estimated parameters as well as in their uncertainty 
estimations. So far we could not decide which explanation is 
more likely. Interestingly this excessive TTV scatter persists 
and is important even in such a robustly fittable case like 
HD 189733, in which the model nonlinearity should be well 
suppressed thanks to a large number of available observa¬ 
tions and lightcurves. 
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Figure 4. Graphs of the TTV residuals, based on the best fitting transit ephemeris (values of P and from Table [4] The transits 
for which a significant red noise was detected in the photometry are labelled with a circle. 
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derived TTV data (306 points) 



normalized residual 


16 points removed 



normalized residual 


Figure 5. Testing the normality of the TTV data derived in the work. The plot is similar to Fig. [3] but refers to the TT V residuals (rela¬ 
tively t o a stric tly periodic transit ephemeris) rather than to the photometric ones. In the right panel, we removed the lSouthworth et al.l 
J2009al lbl 120101) data (7 points), as well as a few potential outliers (9 points). 
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Table 3. Explanation of the TTV data file 


column 

explanation 

i 

Primary data 

integer transit count (number of the transiter’s revo¬ 
lutions, restarts from zero for each star) 

2, 3 

fitted transit mid-time (BJDtdb) and its uncertainty 

4 

Auxiliary data (other parameters of the fit) 
a standardized name of the lightcurve file that includes 
the date, the target name, and the name of the ob¬ 
server or first author of a paper 

5 

the reference time To to which the following trend co¬ 
efficients refer 

6 

number of the following trend coefficients including 
the constant, i.e. the trend degree plus one (always 4 
in this work) 

7, 8 

fitted constant level of the magnitude and its uncer¬ 
tainty 

9, 10 

fitted linear trend coefficient and its uncertainty 
(mag-day -1 ) 

11, 12 

fitted quadratic trend coefficient and its uncertainty 
(mag-day -2 ) 

13, 14 

fitted cubic trend coefficient and its uncertainty 
(mag-day -3 ) 

15, 16 
17, 18 

19 

fitted limb darkening coefficient A and its uncertainty 
fitted limb darkening coefficient B and its uncertainty 
an adopted value of a scale parameter cr sca i e needed to 

20, 21 

fully characterize the regularized model of the photo- 
metric noise, see (|Baluev! 20151) 

fitted photometric white jitter w >,+ fdef. in|Baluevj 

([201^)) and its uncertaintv 

22, 23 

fitted photometric red jitter cr* rec j = y/p re d and its 
uncertainty 

24, 25 

26 

fitted correlation timescale r of the red jitter (days) 
r.m.s. of the best-fit residuals for this lightcurve 












Table 4. Fitted parameters of exoplanetary transit curves 



total 



Assuming fittable TTV 



Fixing TTV residuals at zero 1 


transiter 

number of 

transits 

number of 

red-noised 

transits 

radii ratio 

r = Rpi/R * 

half-duration 

W 2 [days] 

impact par. 

b 

\J~X- TTV 1 

mid-times correl. 

MAD/MAX 

orbital period 

P [days] 

ref. mid-time 

ef [bjdtdb— 

-2450000] 

comment 

WASP-2 b 

38 

9 

0.1355(31) 

0.03701(35) 

0.7380(94) 

1.38 

0.00056/0.0026 

2.15222163(42) 

5894.07919(15) 


WASP-4 b 

43 

20 

0.15495(32) 

0.044898(53) 

0.136(28) 

1.35 

0.00025/0.0039 

1.338231624(68) 

4966.782814(21) 


WASP-5 b 

11 

4 

0.1136(13) 

0.05028(29) 

0.446(40) 

1.34 

0.0070/0.10 

1.62842953(52) 

6446.98868(17) 


WASP-52 b 

22 

7 

0.1629(44) 

0.03858(69) 

0.598(32) 

2.00 

0.0020/0.016 

1.7497835(11) 

6673.82149(13) 


Kelt-1 b 

15 

6 

0.0783(14) 

0.05661(42) 

0.998(44) 2 

2.56 

0.00092/0.017 

1.21749448(80) 

6093.13464(19) 

Showing y/1 — b‘ 2 instead of 
6, see note 2 

Corot-2 b 

20 

7 

0.1639(21) 

0.04736(36) 

1.000(89) 2 

1.17 

0.00079/0.015 

1.74299673(31) 

5628.44758(14) 

Showing \/l — 6 2 instead of 
b , See note 2 

XO-2 b 

25 

10 

0.1036(13) 

0.05563(32) 

0.996(31) 2 

1.56 

0.00056/0.012 

2.61585779(43) 

5139.16092(13) 

Showing y/l — instead of 
6, See note 2 

TrES-1 b 

43 

17 

0.13781(97) 

0.05225(18) 

0.191(67) 

1.54 

0.00038/0.0078 

3.03006973(18) 

5016.969937(70) 

Some of iBakos et alJ 

1200(f) data contain manv 
outliers and possibly 

non-Gaussian noise 

HD189733 b 

67 

29 

0.15712(40) 

0.037576(48) 

0.6636(19) 

2.16 

0.00037/0.071 

2.218575200(77) 

3955.5255511(88) 

GJ436 b 

22 

8 

0.088(10) 

0.02102(71) 

0.806(31) 

1.22 

0.00026/0.0019 

2.64389846(44) 

5280.17568(17) 

Most data are of poor qual¬ 
ity (large trends), and a 
significant eccentricity of ~ 
0.15 is not taken into ac¬ 
count 


^■The orbital periods, the reference mid-times, and the values of Xttv were obtained after removal of the ISouthworth et alJ (|2009al lbl. l20ld l data, as these 
errors. 

2 Impact parameter b is close to zero and is thus a highly nonlinear parameter here. Its estimation is severely nongaussian, and the formal uncertainty cr^ 
handle this pecularity, the values in the column for b are replaced by a = y/1 — b 2 and its uncertainty a a = (b/a)ab, which are more informative here. 


data are affected by clock 
is much larger than 6. To 
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6 RV DATA AND THEIR ANALYSIS 

In addition to the TTV time series from the photometry, 
we also derive RV data for some of our stars using spectra 
found in the HARPS, HARPS-N, and SOPHIE archives. We 
process the se spectra using the advanc ed HARPS-TERRA 
pipeline bv lAnglada-Escude fc Butler! d2012l) . Public spec¬ 
tra were available for the following targets: WASP-2, 4, 5, 
HD 189733, GJ 436, CoRoT-2. 

Most of these RV data appeared less accurate than the 
typical 1 m/s precision demonstrated by HARPS. This is 
because many of these targets are rather faint. Moreover, it 
seems that adequate modelling of these data should involve 
non-trivial treatment of the RV noise. Frequently, these data 
are combined in short series acquired within a few hours or 
even shorter. The data within suc h single-night ser ies should 
be significantly correlated fe.g. lNelson et al.ll2014f ). and this 
effect can be easily detected in some of our data. Usually, 
these short-term runs were clearly intended to catch the 
Rossiter-McLaughlin effect during the planetary transit. On 
larger time scales (days to weeks), the RV noise may still re- 
mai n correlated as well, likely due to the s tellar activity (see 
e.g. lBaluevir2013bl : iRobertson et alJl20l4 l. and this type of 
correlation is different from the one emerging at short time 
scales. 

The full modelling of all these effects is outside of the 
scope of this paper, which was intended to deal mainly with 
the photometry and TTV data. However, these RV data 
may carry useful information that can be helpful in verify¬ 
ing the results of our TTV analysis. In this work we adopted 
a simplified “first look” approach to the RV data analysis. 
First, we replaced the RV series acquired in a single night 
by their averages. This eliminated the need to model the 
Rossiter-McLaughlin effect, as well as the RV noise corre¬ 
lations appearing within a single night. Of course, such a 
procedure likely adds some systematic error to these “cumu¬ 
lative” measurements, but in this work this is a satisfactory 
precision. Then we passed these data throug h a periodog ram 
computing tool of the PlanetPack software (|Baluevll2()13al l. 
taking into account the best fitting contribution from tran¬ 
siting planet (i.e., including it in the null model of the peri¬ 
odogram) . 

In this analysis, none of the periodograms revealed hints 
of any additional variations in the data. One exception is 
the GJ 436 case. In this case we have rather larg e amount o f 
the HARPS data, as well as Keck data from IlKnutson et al.l 
I2014T) . and both these datasets demonstrated clear hints of 
a red noise, similar to the one investigated e.g. by iBaluevI 
(l2013b|j . The same might be true for HD189733, for which 
we also have large amounts of the RV data. We plan to 
consider these cases detailedly in a separate work. 


7 THE WASP-4 CASE 

The periodogram of the WASP-4 TTV data is shown in 
Fig. [fj] and it reveals some marginally significant periodic¬ 
ities, corresponding to the TTV amplitudes ~ 10 — 20 sec. 
However, there is not any stable pattern of the peaks, as 
the periodograms are severely model-dependent. Depending 
on which data subset we include in our anal ysis, we obtain 
different results. Note that the data by iSouthworth et al.l 


d2009bh are not reliable for TTV studies due to the clock 
failures noticed for this telescope. This is why we also con¬ 
sider in Fig.[S]a reduced TTV time series obtained by remov¬ 
ing these data. There are only 4 these TTV points, but they 
have a good formal accuracy as derived from the photomet¬ 
ric fit, and thus affect the TTV periodogram significantly. 
Also, results of the period analysis depend on various other 
subtleties, e.g. on the degree of the polynomial trend used 
in the photometry model. From the other side, although 
these variations are marginal and model-dependent, we can¬ 
not just attribute them to the noise, as we did not observe 
anything similar in the other stars of this study. 

Using the method of Sect. [I] we detecte d an out lier in 
the WASP-4 TTV data, owing to one of the lPetrucci et ahl 
d2013h li ghtcurves. Thi s pecu liar lightcurve was already no¬ 
ticed bv lPetrucci et al.l (1201311 themselves. Interestingly, the 
height of the periodogram peak at ~ 5.14 d is increased after 
removal of this TTV outlier from the analysis. 

The RV data from the HARPS archive did not reveal 
any significant signal in addition to the primary transiting 
planet. However, the star is rather faint for HARPS. Also, 
we should not expect a direct connection between the TTV 
and RV periods, because often TTV is an indirectly induced 
variation, owing to the planetary dynamical perturbations. 
Such a variation may appear due to e.g. a mean motion reso- 
nance, which is not al ways e asily detected in radial velocities 
(e.g. lAnglada-Escude et al.ll201(il f. 

Currently we remain uncertain about the nature of 
these variations that are possibly present in the WASP-4 
TTV data. Possibly, the TTV signal might be more com¬ 
plicated than just a sinusoid, and it should be modelled in 
the framework of the Newtonian N-body fits. Another possi¬ 
bility is star spots inducing systematic perturbations in the 
timing measurements. In any case, it seems that we should 
keep tracking this target or maybe even focus increased at¬ 
tention on it by making further observations in future. 


8 CONCLUSIONS AND DISCUSSION 

As the main topic of the manuscript is the TTV exoplanet 
detection, it is interesting to investigate its efficiency rela¬ 
tively to the more classic methods like the Doppler planet 
detection. An apparent TTV signal can be induced via two 
mechanisms: the dynamical perturbation on the transiter’s 
motion and the light arrival time delay due to the finite light 
speed (the Roemer effect). The dynamical perturbations are 
difficult to predict in a general case, as they severely depend 
on many orbital parameters (e.g. they may drastically in¬ 
crease when planet move in a mean-motion resonance). But 
the Roemer delay can be easily assessed, so let us now con¬ 
sider it detailedly. 

Assume that a distant second planet has the mass of m 
and orbits the star (with the transiter) on a circular orbit 
with the radius of a. Such a planet should induce a similar 
circular motion on the host star on an orbit of the radius 
a' ~ m/M+a, where A/* is the star mass. This reduces to a 
sinusoidal variation of the transits time delay with an am¬ 
plitude of 

Attv — — T-r- sini, (11) 

cAZ* 
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period [day] period [day] period [day] 

Figure 6. Periodograms of the WASP-4 TTV data, showing a few marginally signi ficant peaks. The plot i n the left panel was based on 
all available transit lightcurves, the middle panel is for a reduced data set with the lSouthworth et al.H2009bh data removed, and in the 
right panel one TTV outlier detected in the jPetrucci et al.|[2013h data was also removed. The thick horizontal lines label the significance 
levels of lcr or FAP « 31.7 per cent, 2<r or FAP « 4.6 per cent, 3rr or FAP « 0.27 per cent. The significance and the period for the tallest 
peak are printed in each panel. All periodograms involve a quadratic trend of the TTV delay with fittable trend coefficients. 


where c is the speed of light, and i is the orbital inclination 
to the sky plane. 

In the same case, the amplitude of the Doppler variation 
induced on the star is equal to 


Kbv — m sin i 



( 12 ) 


where G is the gravitational constant. The ratio of these 
amplitudes now looks like 

7 3/2 

(13) 
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where P is the orbital period of the distant companion, com¬ 
puted according the the third Kepler law. We can see that 
the relative efficiency of the TTV method, in comparison 
with the RV one, only depends on the orbital period of the 
unseen companion. 

Contrary to the RV and classic transit surveys, the TTV 
method might be useful to detect distant companions. Of 
course, in any case we must track the variation over at least 
a single period, so detection of a distant companion neces¬ 
sarily requires a long observation run. From this point of 
view, ground-based observatories should be more useful for 
TTV planet detection. Spacecraft rarely operate over a term 
longer than a few years, while ground-based observations can 
run on an indefinitely long time base. 

The accuracy of the TTV data presented in this work 
is such that it would allow a robust detection of a signal, 
if its amplitude is above 1 min. With this TTV threshold, 
the formula m implies that a Jupiter-mass planet orbiting 
a Solar-mass star could be detected, if its semimajor axis 
was at least 60 AU, implying an extremely large period of 
at least ~ 400 yr. This is not a realistic requirement for the 
observation time. Alternatively, the planet mass should be 
at least 12 times the Jupiter mass, if we want to detect it 
at a Jupiter-like orbit with a = 5.2 AU. Also, the exoplanet 
detectability naturally increases for small-mass hosts. Thus, 
to robustly detect a Jupiter-mass planet, we should either 
wait for a long time or to decrease the TTV measurement 
errors to a level of seconds, roughly by an order of mag¬ 
nitude. Based on these computations, we may say that the 
TTV detection method is now in its early development stage, 
comparable to the early era of the Doppler technique (prior 
or near the 51 Pegasi b detection). 

Summarizing, we can conclude that amateur-class TTV 
observations may occupy two possible niches: (i) detection 


of long-period massive exoplanets and (ii) detection of ex¬ 
oplanets trapped in mean motion resonances with known 
transiters (thus generating a more remarkable TTV signal 
via dynamical effects). However, the TTV detection of long- 
period Jupiter-like planets might be a feasible task for more 
advanced ground-based observatories that can achieve the 
TTV accuracy of a few seconds. Some of the data listed in 
Table[l]do provide such an accuracy. To enable the detection 
of Jupiter twins, such observations should be carried out in 
an experimental monitoring regime, i.e. on a regular basis 
and over a long term, similar to the modern Doppler sur¬ 
veys. Note that Jupiter analogues may represent a special 
interest, because such planets would point out exoplanetary 
systems that have architectures similar to Solar System, in 
which all giant planets are quite distant from the Sun. The 
chance to find an Earth twin in such a system might be 
higher. 

At last, we note that the transit fitting algorithm pre- 
sent ed here is now i mplemented in the free PlanetPack pack¬ 
age dBaIuevH2013aT ) and is made available as of the current 
version PlanetPack 2.1. 
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